Probing the exchange field of a quantum-dot spin valve by a superconducting lead 
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Electrons in a quantum-dot spin valve, consisting of a single-level quantum dot coupled to two 
ferromagnetic leads with magnetizations pointing in arbitrary directions, experience an exchange 
field that is induced on the dot by the interplay of Coulomb interaction and quantum fluctuations. 
We show that a third, superconducting lead with large superconducting gap attached to the dot 
probes this exchange field very sensitively. In particular, we find striking signatures of the exchange 
field in the symmetric component of the supercurrent with respect to the bias voltage applied 
between the ferromagnets already for small values of the ferromagnets' spin polarization. 
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I. INTRODUCTION 

The rapidly evolving field of spintronics^'^ not only 
takes into account the charge degree of freedom of elec- 
trons as does conventional electronics, but additionally 
makes use of the spin degree of freedom thereby open- 
ing up interesting new applications as well as address- 
ing questions of fundamental research. Quantum dots 
coupled to ferromagnetic electrodes provide one partic- 
ular example of spintronic systems. They can be ex- 
perimentally realized in a number of ways, e.g., us- 
ing metallic nanoparticles,'^^^ semiconductor quantum 
dots,'"^^ quantum dots defined in nanowires,^^ carbon 
nanotubes,^'^"^'' or single molecules. 

From the theoretical point of view, many different 
properties of such systems have been studied, e.g., spin- 
diode behavior in quantum dots coupled to one ferro- 
magnetic and one normal lead,^^'^° the angular depen- 
dence of conductance, ^^"^^ negative tunnel magnetoresis- 
tance,^^"'^^ and current noise. ^'^"'^^ Of particular interest 
are quantum dots coupled to noncoUinearly magnetized 
ferromagnets. These systems on the one hand show a 
nonequilibrium spin accumulation on the quantum dot 
that has the tendency to block transport. On the other 
hand, there is an effective exchange field acting on the 
dot spin that is caused by virtual tunneling between the 
quantum dot and the leads and gives rise to a precession 
of the spin accumulated on the dot.^^'^^ The interplay 
between these two effects gives rise to a number of in- 
teresting transport properties as, e.g., a deviation from 
the harmonic dependence of the linear conductance on 
the angle enclosed by the magnetizations,^^ a u-shaped 
conductance curve with a broad region of negative dif- 
ferential conductance,^^ a nontrivial bias dependence of 
the Fano factor and characteristic features in the finite- 
frequency noise at the Larmor frequency associated with 
the exchange field'^^ as well as to a splitting of the Kondo 
resonance. '^^"'^^ Detecting the exchange field in experi- 
ments is quite challenging as most of the effects listed 
above rely on a strong spin blockade of the quantum dot 
that exists only for large polarizations of the leads. For 



this reason, the exchange field has so far been detected 
only by the induced splitting of the Kondo resonance 
in Cgo molecules, InAs quantum dots,^'^^ and carbon 
nanotubes^*^ coupled to ferromagnetic leads, respectively. 

Here, we propose an alternative way to experimen- 
tally access the infiuence of the exchange fields on the 
transport properties in the regime of weak tunnel cou- 
pling by adding a superconducting electrode to the quan- 
tum dot as shown schematically in Fig. 1. Quantum 
dots coupled to superconducting electrodes are interest- 
ing on their own as they show an interplay between strong 
Coulomb interaction that has the tendency to destroy 
superconducting correlations on the dot and nonequi- 
librium effects that can help to induce a superconduct- 
ing proximity effect on the dot. In the subgap regime, 
transport between the quantum dot and superconduc- 
tor takes place via Andreev refiections which have been 
investigated theoretically extensively^"*"^^ Further stud- 
ies involved multiple Andreev reflections^'^'^* and trans- 
port in the Kondo regime. ^^^^^ Experimentally, quan- 
tum dots coupled to superconductors have been realized 
using carbon nanotubes,^^"^* graphene,^^ semiconductor 
nanowires,^^'^*^'^^ self-assembled semiconductor quantum 
dots,^^ and single molecules.'''^ 

For the system under investigation, we compute the 
current into the superconductor using a real-time dia- 
grammatic approach^^'^""^^'^"'"^^ in the limit of an infi- 
nite superconducting gap. This is a reasonable approxi- 
mation as long as the excitation energies of the quantum 
dot are smaller than the gap such that sub-gap transport 
takes place. The current has, in general, even and odd 
components with respect to the voltage applied between 
the ferromagnets. We find that the even component ex- 
hibits clear evidences of the exchange field even for small 
polarizations of the ferromagnets. 

The paper is organized as follows. We introduce the 
model for the quantum dot coupled to the three elec- 
trodes in Sec. II. In Sec. Ill, we describe how to apply 
the real-time diagrammatic theory to the system under 
investigation. We present our results in Sec. IV and give 
a conclusion in Sec. V. 
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FIG. 1. (Color online) Schematic model of a quantum-dot 
spin valve with an additional superconducting electrode at- 
tached. A single-level quantum dot with excitation energies 
e and e + U is coupled to to noncoUinearly magnetized ferro- 
magnets and a superconductor via tunnel barriers. 



II. MODEL 

We consider a quantum-dot spin valve, i.e., a quantum 
dot coupled to two ferromagnetic electrodes with mag- 
netizations pointing in arbitrary directions riL, hr. In 
addition to the two ferromagnetic electrodes we consider 
a third superconducting lead tunnel coupled to the quan- 
tum dot. The Hamiltonian of the system hence consists 
of different terms describing the ferromagnetic electrodes, 
the superconducting electrode, the quantum dot and the 
tunnel coupling between the dot and the leads, 



H 



Hr + Hs + -ffdot + -ffti 



(2.1) 



We describe the ferromagnetic electrodes with chemical 
potentials fir, r = L, R, by means of the free-electron 
Hamiltonians 



Hr 



^(e-rkcr 



(2.2) 



where the quantization axis is chosen to be parallel 
to the magnetization of the respective lead. We as- 
sume the densities of states Prai^^) — J^u^i'^ ~ ^rka) 
to be constant, prcr(<^) — Pra, and spin dependent. 
The asymmetry between majority {a — +) and minor- 
ity (cr = — ) spins is characterized by the polarization 
Pr — iPr+ — Pr-) / iPr+ + Pr-) which varics between p = 
for a nonmagnetic electrode and p = I for a halfmetallic 
electrode with majority spins only. 

We model the superconductor by means of a mean-field 
BCS Hamiltonian with a superconducting gap A, which 
can be chosen to be real and positive without loss of gen- 
erality. We choose the chemical potential of the super- 
conductor as reference for energies and set it to zero. In 



the limit of an infinite superconducting gap, A — oo, the 
single-level quantum dot tunnel coupled to the supercon- 
ductor is described by the effective dot Hamiltonian''^"^'' 



dot.eff 



ec. 



A t t 
9 "^t^i 



where e is the energy of the spin-degenerate level in the 
dot and U denotes the Coulomb energy for double oc- 
cupancy of the dot. The effective pair potential Fg in 
Eq. (2.3) is the tunnel-coupling strength between the dot 
and the superconductor and it is related to microscopic 
parameters by Fg = 27r|tgppg, where ig is the tunnel- 
ing amplitude between dot and superconductor and pg is 
the normal-state density of states of the superconducting 
lead. 

The effective dot Hamiltonian accounts for the cou- 
pling to the superconductor exactly. This allows us to 
deal with an arbitrarily strong superconductor-dot cou- 
pling Fg. The eigenstates of the effective dot Hamilto- 
nian (2.3) are given by the singly occupied states | f) and 
I \) as well as by the two states \+) and |— ). The latter 
ones are linear combinations of the empty and doubly 
occupied dot states |0) and \d) = clc||0): 




(2.4) 



(2.5) 



The energies of the eigenstates are given hy ^ ~ e 
and E± = 5/2 ± ea. Here, 5 = 2£ -f C/ denotes the 
detuning from the particle-hole symmetry point while 
2eA = Y^(5^ -I- F| measures the energy difference between 
the states |-|-) and |— ). 

We define the Andreev bound-state energies as the ex- 
citation energies of the dot in the absence of tunnel cou- 
pling to the ferromagnets: 



7 



7 



F2 

S' 



7,7 



(2.6) 



2 2 

As we allow arbitrarily oriented magnetizations of the 
ferromagnetic leads, riL, nR , it turns out to be most con- 
venient to quantize the spin on the quantum dot in the 
direction of iil x n^. In this case, the tunnel coupling 
between the dot and the ferromagnets is characterized by 



tun,F 



E 



+«rk- (-e'^'-ct + e-'^-c;)l + h.c, (2.7) 



i.e., the majority/minority spin electrons of the leads cou- 
ple to both spin up and spin down states of the quantum 
dot. In the tunnel Hamiltonian, 0l — ~(t>K — (t>/2 is 
half the angle between the magnetizations. The tunnel 
matrix elements tr can be related to the spin-dependent 
tunnel couplings 27r jf^P Pro- ■ The total tunnel coupling 
to lead r is given by F^ = J2a '^'^\tr\'^Prcr/'^- 
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III. TECHNIQUE 

In order to compute the transport properties of the 
quantum-dot spin valve with an additional supercon- 
ducting lead, we use the diagrammatic real-time tech- 
nique^*"''^ in its extension to ferromagnetic^^ and super- 
conducting^^"^^ leads. 

The basic idea of this approach is to integrate out the 
noninteracting, fermionic degrees of freedom of the elec- 
trodes. The remaining system consisting of the quantum 
dot is then described using a reduced density matrix p'^°'^ 
with matrix elements P^^ = (X2|p"'^lxi)- For conve- 
nience, we introduce for the diagonal matrix elements 



the abbreviation 



px 

X 



For a quantum dot coupled to both ferromagnetic and 
superconducting leads the reduced density matrix in the 
stationary limit takes the form 



P+ \ 
P_ ' 




(3.1) 



V p| p^y 



in the basis |-|-), |— ), | t), and | I). We note that while 
coherent superpositions of states | t) and | l) have to be 
taken into account as we are dealing with a noncoUinear 
geometry, coherent superpositions of |-|-) and |— ) can be 
neglected if Fs ^ F^ (which we assume from now on) as 
they are at least of order 0{Tr ) and therefore do not con- 
tribute to transport in the limit of weak tunnel coupling 
between dot and ferromagnets. 

In the stationary state, the reduced density matrix 
obeys a master equation of the form 

= P^^ = -^(i?,, - E^,)P^f + J2 W^Pi- (3-2) 

X'lX2 

While the first term on the right-hand side describes the 
coherent evolution of the system, the second term char- 
acterizes the dissipative coupling to the electrodes. The 
quantities W^^^? are generalized transition rates in Li- 
ouville space. For tunnel couplings small compared to 
the temperature, F,. <C k^T, they can be evaluated in a 
perturbation expansion in the tunnel couplings as irre- 
ducible self-energy blocks of the dot propagator on the 
Keldysh contour. The corresponding diagrammatic rules 
are summarized in Appendix A. 

Using the effective dot description and introducing the 
average spin on the dot as 



P\ 



P 



1 Sy 



pi 



P 



2i 



(3.3) 



as well as the probability Pi = P-^ + P^ to find the 
dot singly occupied, we can split the master equation 
into one set for the occupation probabilities and one set 
for the average spin. We introduce the abbreviations: 

r.± = F, (li^l^) and = /±(£;a,77'), where 



/+ (w) = 1 — f~ (w) denotes the Fermi function of lead 
r. The equations governing the occupation probabilities 
are given by 



= 



dt 



p =Ea. p 



Pi 



Pi 



Prhr{Sr ■ n^), 

(3.4) 

where the expressions for the matrices A,, and the vec- 
tors hr are given in Appendix B. We see that similar to 
the case of an ordinary quantum-dot spin valve, the 
dynamics of the occupation probabilities couples to the 
average spin accumulated on the quantum dot. 

As in the case of the normal quantum-dot spin valve, 
the master equation for the average dot spin can be cast 
into the form of a Bloch equation. 



dS _ /dS\ 



where the first term 



dS\ 



dS 

Itt 



(3.5) 



dS\ 
'dt „ 



2 ^ 



[(-F,+/-+++F,_/+__)P+ 



+ (-r.-/-+_ + r,+/++)p_ 

+F,+ /,+ + + F,_/++_) Pi] PrUr, (3.6) 

describes the nonequilibrium spin accumulation on the 
dot due to spin-dependent tunneling of electrons onto 
the dot. The relaxation of the dot spin is described by 
the second term. 



dt I 



2 ^ 



-F,+ /++++F,_/,V)S, (3.7) 



which is proportional to the spin accumulated on the dot. 
The dot spin relaxes either by electrons with a given spin 
leaving the dot to the ferromagnetic leads or by electrons 
with a spin opposite to that on the dot entering the dot 
from the ferromagnets, thus forming a spin singlet. Fi- 
nally, the last term 



dS 



= X B^, 



(3.8) 



prcc 



describes a precession of the dot spin due to an exchange 
field which is given by 



Br 



Pr^r 

27r 



^ 7'F^^ Re * 



1 . /3(£'a.7'7 - Hr) 
2+* 277 



(3.9) 

where 4* is the digamma function. The exchange field 
is the manifestation of a spin-dependent level renormal- 
ization due to virtual tunneling between the dot and the 
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FIG. 2. (Color online) Amplitude of the exchange field Br-n,. 
in units of PrTr for k^T — O.Olf/, Fs = OAU as a function 
of the chemical potential Hr and detuning 5. The peaks and 
dips map out the Andreev bound states whose energies are 
indicated by dotted lines. 



ferromagnetic electrodes. We emphasize that the cou- 
phng to the superconductor influences the exchange field 
only through the position of the Andreev bound states. 
As can be seen in Fig. 2, the exchange field takes on 
large values whenever one of the Andreev bound states 
is in resonance with the Fermi level of the ferromagnet. 
This behavior is similar to the ordinary quantum-dot spin 
valve where the exchange field becomes maximal at res- 
onance as well. 

The particle current flowing from the ferromagnetic 
leads into the quantum dot is given by 

ir=Y: w'^iilP^.- (3-10) 

XX1X2 

Here, W^''^^^ are the current kernels that can be ob- 
tained from the generalized transition rates W^^^ by 
multiplying the rate with the number of electrons trans- 
ferred from lead r to the quantum dot in the process 
associated with this rate. We flnd 

Ir = (r,r-/,t — - rr+fr+ + ) P+ 

+ (r.+./+ + -r,_/,-_)p_ 

+ 2 (^r+fr++ + ^r-fr+- 

-r.-/-__-r,+/-_+)Pi 

~Pr {Tr-fr— + ^r+fr- + 

+Tr+f+++ + r,_/++_) S • n,. (3.11) 

In the stationary state, the current into the supercon- 
ductor is related to the currents between the dot and the 
ferromagnets by current conservation, which is automat- 
ically satisifed in the real-time diagrammatic theory, 

/s=/l + /r. (3.12) 
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FIG. 3. (Color online) Transport processes in the proximized 
quantum-dot spin valve. In the first process, the |— ) state 
is projected onto the |0) state, pushing a Cooper pair into 
the superconductor as indicated by the double arrow. Then 
an electron tunnels in from the left lead, leaving the dot in 
the singly occupied state. Similarly, in the second process, 
an electron leaves the singly occupied dot to the right lead. 
Then, a Cooper pair enters from the superconductor to bring 
the dot in the state |— ). As both processes probe the empty 
contribution to |— ), their rates are proportional to l-|-5/(2eA). 
Similarly, the other two processes probe the doubly occupied 
component of |— ) such that their rates are proportional to 
1-(5/(2£a). 



Experimentally, one can therefore measure either the cur- 
rent flowing into the superconductor or the difference be- 
tween the currents that enter from the left and leave to 
the right ferromagnet. In the following discussion we will 
however always refer to the current into the superconduc- 
tor for simplicity. 



IV. RESULTS 

In this section, we discuss how the current into the 
superconductor can be used to probe the exchange field. 
We consider symmetrically biased ferromagnets, ~ 
— /iR = fjL, while the superconductor is kept at = 0. 
We split the current into a component that is a symmet- 
ric function of bias, /g-*"°(/x) — [/s(m) + -fs(— m)]/2, and 
a component that is an antisymmetric function of bias, 
jantisym^^^ _ j/g^^") _ /g(— /i)]/2. Thc Symmetric com- 
ponent of the current turns out to be very sensitive to 
the exchange field. 



A. Symmetric quantum-dot spin valve 

We start our discussion by considering a symmetric 
system, i.e., we assume the tunnel couplings to the left 



5 



and right ferromagnet to be equal, Fl = Fr = F/2. Fur- 
thermore, we assume both ferromagnets to have the same 
polarization p. In this particular case, due to symmetry, 
the current possesses only a symmetric component with 
respect to /i, i.e. Is{^J■) = Is''^{^J■)■ 

We first explain why the supercurrent vanishes for 
coUinear geometries. We then show that the spin ac- 
cumulation in the noncoUinear configuration gives rise 
to a finite current into the superconductor that is sensi- 
tive to the exchange field. Finally, we show that a spin 
relaxation on the dot reduces the supercurrent but nev- 
ertheless still allows a detection of the exchange field. 

In general, the current into the superconductor van- 
ishes in the small-bias regime Ea h < M < , where 

the quantum dot is Coulomb-blockaded. For a symmetric 
system, the current also vanishes in the large-bias regime, 

fi > Ea++ or /i < Ea , where all dot states contribute 

to transport, due to particle- hole symmetry. Hence, we 
can expect a finite current into the superconductor only 

in the intermediate bias regime, Ea < fi < Ea y or 

Ea+- < fJ, < Ea++- According to Eq. (3.11) and (3.12), 
the supercurrent in this regime is given by 



0.001 



/s = 



-2 — P+ 
£A 



— p_ Pi 

Ea 2eA 



'pS ■ (riL + hr) 



pS 
2^ 



S • (riL - hr) 



(4.1) 



For parallel magnetizations, where the spin accumulation 
on the dot vanishes and the dot occupation probabili- 
ties satisfy Pi = 2P_ = 1/3, we find that the supercur- 
rent vanishes also in the intermediate bias regime. To 
understand the mechanism behind this behavior, let us 
consider the transport processes that contribute to the 
superconductor. Notice that in the intermediate regime 
the state |-|-) is inaccessible. We find that the contribu- 
tions from the first two processes shown in Fig. 3 cancel 
each other. They both transfer equal amounts of charge 
between the dot and the superconductor when project- 
ing the state |— ) onto the state |0). Furthermore, both 
processes have identical rates (a factor of 2 due to spin 
is compensated by Pi — 2P_). In consequence, they give 
rise to a vanishing supercurrent for any value of the de- 
tuning S. Similarly, one can show that the other two 
processes that probe the doubly occupied component of 
|— ) cancel each other. 

The situation is more complex in the antiparallel con- 
figuration due to the finite spin accumulation on the 
quantum dot. Processes 1 and 4 in Fig. 3 which build 
up the spin accumulation have rates proportional to 
1 + (5/(2£a) and 1 — 5/{2eA), respectively. Hence, in sum 
the spin accumulation is insensitive to the detuning S. A 
similar reasoning holds for processes 2 and 3 that relax 
the dot spin. As the supercurrent vanishes a.t S — due 
to particle-hole symmetry, it therefore has to vanish for 
all values of 6. 

For noncoUinear geometries, we find a finite current 
into the superconductor, cf. Fig. 4. To understand the 
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FIG. 4. (Color online) Current into the superconductor for a 
symmetric coupling, Fl ~ Fr,, /il = —MR = M perpen- 
dicularly magnetized ferromagnets, = 7r/2. In the upper 
panel, we have p = 0.3, while in the lower panel p — 0.95. 
The solid (black) curves take into account the exchange field, 
it is neglected in the dashed (red) curves. Other parameters 
are S = 0.2U, Fs = 0.4(7, ksT = O.OIC/. 



mechanism leading to the finite current, we first neglect 
the exchange field in our discussion and turn to its effect 
afterwards. 

In contrast to the antiparallel configuration, in the 
noncoUinear geometries the spin accumulation on the 
quantum dot (and therefore also the probability to find 
the dot singly occupied) is sensitive to 6. This can be 
understood by considering again the processes 1 and 4 of 
Fig. 3 which are responsible for the spin accumulation. 
While process 1 builds up a dot spin in in the direction of 
+nL, process 4 build up a spin in the direction of — hr . 
For 6 = both processes contribute equally to the spin 
accumulation. We therefore find that a spin builds up 
in the direction of riL — hr. For positive 5, process 1 
dominates and hence the spin points toward riL, while 
for negative S, process 4 dominates and the spin points 
toward — hr. As for finite S the left-right symmetry is 
broken, the cancellation between the supercurrent con- 
tributions from processes 1 and 2 (3 and 4, respectively) 
does not hold any longer and a finite supercurrent can 
fiow. Neglecting the exchange field, the following ana- 
lytic expression for the current into the superconductor 
can be found: 



/s = 



EA (48ei - 2p2(2,52(l + 2cos0) - F2(1 - cos </))))■ 



(4.2) 

The above formula shows that the current into the super- 
current in the noncoUinear geometry flows for any value 
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FIG. 5. (Color online) When the third lead is a normal metal 
instead of a superconductor, a finite current flows for any 
magnetic configuration, thereby completely obscuring the ex- 
change field effects for small polarizations. Parameters are 
5 = OAU, Fr = Tn, p = 0.3, ksT = O.OIU. 
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FIG. 6. (Color online) Influence of spin relaxation on the 
supercurrent in the noncoUinear configuration. As the spin 
accumulation is reduced, the supercurrent is decreased. How- 
ever, the exchange field effects still prevail. Polarization is 
p = 0.95, other parameters as in Fig. 3. 



of the polarization p ^ 0. It is only the magnitude of the 
current that is affected by the strength of the polariza- 
tion. Using realistic parameters, p = 0.3, U ~ ImeV, 
Fs ^ 0.5 meV, F ^ 100/LteV, and 4> = 7i'/2, we obtain as 
an order of magnitude of the current /s ~ 1 pA which is 
challenging but not impossible to measure with current 
techniques. 

If the exchange field is taken into account, there is still 
a finite current flowing in the intermediate bias regime. 
As the dot spin now precesses in the energy-dependent 
exchange field, it acquires a finite z component while si- 
multaneously the X and y components which influence 
the supercurrent, cf. Eq. (4.1), get reduced and show 
a nontrivial bias dependence. In consequence, the su- 
percurrent also deviates from its steplike behavior in the 
absence of the exchange field and even changes sign. Fur- 
thermore, there is a finite supercurrent flowing in the 
large bias regime because the symmetry-breaking spin 
accumulation on the dot persists in this regime. 

The nontrivial bias dependence of the supercurrent 
opens up the possibility to detect the exchange fleld ex- 
perimentally, even for small polarizations. This is in 
strong contrast to the other exchange field effects that 
arise in the sequential tunneling regime, as, e.g., negative 
differential conductance^^'^^ or the nontrivial dependence 
of current, Fano factor and higher current cumulants on 
the angle between the magnetizations. ^^■•^^'^^ While for 
quantum dots that couple only to ferromagnetic leads all 
exchange field effects rely on a strong spin blockade, in 
the system under investigation here, it is the cancella- 
tion between different transport processes involving the 
superconductor that provides the necessary sensitivity to 
the spin accumulation and the exchange field. To illus- 
trate this, let us consider a system where the supercon- 
ductor is replaced by a normal metal that is coupled to 
the quantum dot with coupling strength Fn- In Fig. 5 
we show the current into the normal metal evaluated to 



first order in F^, an approximation valid if Fn ^ k^T. 
We find that now there is indeed a finite current for all 
magnetic configurations because for a positive (negative) 
detuning only processes where electrons leave (enter) the 
dot to (from) the third lead are possible. In consequence, 
the exchange field effects become practically invisible for 
small polarizations as they are obfuscated by the large 
background current. 

Finally, we discuss the effect of an intrinsic spin re- 
laxation on the dot which we model by adding a term 
— S/r to the right-hand side of the spin master equation, 
Eq. 3.5. Possible mechanisms for such a spin relaxation 
are the coupling to nuclear spins in the quantum dot*^"^^ 
or spin-orbit interaction on the dot.®^'^^ In Fig. 6, we 
show the current into the superconductor for different 
values of the relaxation time r. As the relaxation time is 
decreased, the current is reduced in agreement with our 
discussion above which showed that the spin accumula- 
tion on the dot is crucial to get a finite current. However, 
we also notice that the exchange field effects still remain 
visible when considering a finite relaxation. This shows 
that an experimental detection of these effects should be 
feasible. 



B. Asymmetry effects 

We now turn to the discussion of the situation where 
Fl Fr. We parametrize the tunnel couplings as 
Fl = (1 + a)F/2 and Fl = (1 - a)F/2 such that the 
parameter a with —1 < a < +1 characterizes the degree 
of asymmetry. In this case, the antisymmetric compo- 
nent of the supercurrent with respect to the applied bias, 
ig (/i), IS m general non vanishing. 

We find that for a ^ a finite supercurrent arises in the 
intermediate and large bias regime for all magnetic con- 
figurations. For parallel and antiparallel magnetizations 
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in the intermediate bias regime, while it is given by 
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FIG. 7. (Color online) Symmetrized and antisymmetrized 
current into the superconductor as a function of bias voltage 
for different magnetic configurations and small {p — 0.3, up- 
per panel) and large {p — 0.95, lower panel) polarizations. 
The asymmetry is a = 0.05, other parameters as in Fig. 4. 
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(i-p2)52 + r|' 
(i-aV)^2 + r| 



(4.5) 
(4.6) 



in the large bias regime. In these formulas we assumed 
the current to flow from the left to the right. For a 
current in the opposite direction, one has to substitute 
a — >■ —a. Since the corresponding formulas for the non- 
coUinear case are rather lengthy, we do not give them 
here. 

From the above formulas we read off that the total su- 
percurrent increases as the asymmetry is increased. Fur- 
thermore, we find that the current is decreased when the 
polarization is increased. This means that for experimen- 
tally realistic polarizations an asymmetry can give rise to 
a background current that dominates over the exchange- 
field signal. This problem can be overcome by looking 
at the symmetric and antisymmetric components of the 
current with respect to the voltage. 

In Fig. 7, we show these two quantities as a function 
of applied bias for small and large polarizations. For 
small polarizations, we find that the antisymmetric con- 
tribution shows a steplike behavior that does not reveal 
any exchange field effects and is nearly insensitive to the 
magnetic configuration. In contrast, the symmetric con- 
tribution again reveals the characteristic peaks and dips 
that we encountered already in the symmetric system and 
that are a clear indication of the exchange field. For large 
polarizations, the system behaves rather similarly. The 
symmetric current contribution shows clear signs of the 
exchange field while the antisymmetric part is dominated 
by current steps. However, we now find that also the 
antisymmetric contribution is sensitive to the exchange 
field. 

Hence, we have seen that an asymmetric coupling to 
the ferromagnets gives rise to a finite supercurrent for 
all magnetic configurations that could make the experi- 
mental detection of the exchange field effects difficult for 
small polarizations. To overcome this obstacle, we pro- 
pose to measure the current symmetrized with respect to 
the bias voltage as this allows to recover the exchange 
field effects. 



V. CONCLUSIONS 

magnetizations, respectively, the supercurrent is given by i j ^ ^ u ^ j i ■ 

° J f jj f a J YVe analyzed transport through a quantum-dot spm 

valve with an additional superconducting electrode. We 
find that in the case of noncoUinear magnetization even 
p 2aTg{l — p'^){2eA — aS) for small polarizations, the symmetric component of the 

s ~ £a [(3 + a^P^)rg — (1 — p'^)6{aeA ~ (3 ~ a^)'^)] supercurrent with respect to the applied bias voltage ex- 

(4 3) hibits strong exchange-field effects. In particular, for a 
2/ 2^ system that couples symmetrically to the ferromagnets 

/AP = p " si P_2 (4-4) which are at opposite bias, the supercurrent has only a 

£a [(1 + 3p2)ai5 -|- 2(3 4- p2)eA] symmetric component in bias voltage. In this case, a 



8 



finite supercurrent can only flow for noncollinear magne- 
tizations, as this configuration breaks the left-right sym- 
metry for finite detuning 6. Due to the presence of an 
exchange field acting on the dot spin in noncollinear ge- 
ometries, the supercurrent exhibits a nontrivial bias de- 
pendence and even changes sign. Interestingly, these ef- 
fects occur for any polarization of the ferromagnets p > 0. 
Furthermore, they are robust towards a relaxation of the 
dot spin. For systems with different couplings to the 
ferromagnets, the supercurrent becomes finite for any 
magnetic configuration. We find that for small polar- 
izations, the contribution due to the asymmetry of the 
system dominates over the one due to noncoUinearity. 
We show, however, that by considering the supercurrent 
symmetrized with respect to the applied bias voltage, one 
can extract the exchange field effects also in this case. We 
therefore propose to experimentally access the exchange 
field by measuring the bias dependence of the supercur- 
rent. 
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Appendix A: Diagrammatic rules 

In this appendix, we give the diagrammatic rules (see 
also Ref. 22, 50-52, 74-77) necessary to compute the ker- 
nels W^^^'^ and W^''^^^'^ that occur in Eq. (3.2) and 

XlXi XlXi ^ ^ ' 

Eq. (3.10) in the limit of an infinite superconducting gap, 
A oo. 



3. The tunneling lines involving the ferromagnetic 
electrode r give rise to a factor of ^f^{uji) if they 
do not flip the spin of the tunneling electron. If 
they flip it from up to down, they give rise to a 
factor of l^e"^"^/^ (wi). Flipping the spin in the 
opposite direction gives rise to the complex conju- 
gate of the aforementioned factor. Here, the upper 
(lower) sign refers to lines running forward (back- 
ward) with respect to the Keldysh contour. 

4. Associate with each vertex that annihilates (cre- 
ates) a dot electron with spin a a factor (X2|ccr|xi) 
((X2|cJ.|xi))- Here Xi and X2 are the states that 
enter and leave the vertex, respectively. 

5. Assign an overall prefactor (— z)(— 1)'^+'' where a is 
the number of vertices on the lower propagator and 
b is the number of crossings of tunneling lines. 

6. Integrate over the energies of the tunneling lines 
and sum over all diagrams. 

7. To obtain the generalized current rates W 

multiply each rate W^'^^? with a factor for all tun- 
neling lines that are associated with lead r that is 
the sum of the following numbers: 

(a) -|-1 if the line is going from the lower to the 
upper propagator, 

(b) —1 if the line is going from the upper to the 
lower propagator, 

(c) else. 

Appendix B: Equation for the occupation 
probabilities 



1. Draw all topological different diagrams with ver- 
tices on the propagators. Assign states x and cor- 
responding energies to the corners and all prop- 
agators. The vertices are contracted pairwise by 
tunneling lines that either conserve or flip the spin. 

2. Assign to all diagrams a resolvent l/(A£^ + zO+) for 
each section on the contour between two adjacent 
vertices. Here IS.E is the energy difference between 
the left- and right-going propagators and tunneling 
lines. 



In this Appendix we give the expressions for the quan- 
tities appearing in Eq. (3.4) which can be obtained as 

linear combinations of the kernels W^'^^? by reformulat- 

XlXi 

ing the master equation (3.2) in terms of the occupation 
probabilities and average dot spin. We employ the abbre- 
viations: r^± = F,. {\ ± and f^^^, = f^{EA^jy), 

where f^i^) = 1 ^ fri'^) denotes the Fermi function of 
lead r. 

The matrix depends on the position of the Andreev 
excitation energies and it reads 



^r+fr++ ^r-f^- — j fr + 2^ /r++ 



^r+fr++ + Tr-Zj, ^r-Ir-\ — + ^r+fr — h 2~/r 2~ 'A- — h 2~ ■fr++ 2~ -frA — 



I 

The vectors describing the influence of the average spin on the dot on the diagonal probabilities are given by 

^r+f, — I- ^ r,._/,,^ — 

^^r—fr ^ ^r+fr — \- + ^r+,fr++ + — .fr+- 

(B2) 
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